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We perform a critical assessment of the accuracy of DFT-based methods in predicting stable 
phases within the Co-Pt binary alloy. Statistical mechanical analysis applied to zero kelvin DFT 
predictions yields hnite-temperature results that can be directly compared with experimental mea¬ 
surements. The predicted temperature-composition phase diagram is qualitatively incompatible 
with experimental observations, indicating that the predicted stability of long-period superstruc¬ 
tures as ground states in the Co-Pt binary is incorrect. We also show that recently suggested 
methods to better align DFT and experiment via the hybrid functional HSE06 are unable to resolve 
the discrepancies in this system. Our results indicate a need for better verihcation of DFT based 
phase stability predictions, and highlight fundamental flaws in the ability of DFT to treat late 3d-5d 
binary alloys. 
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I. INTRODUCTION 

Density functional theory (DFT) has grown to be¬ 
come the most popular electronic structure calculation 
method to datei. Modern computational resources have 
made DFT viable as a high-throughput materials design 
technique^Ti^, whereby the existence, stability, and prop¬ 
erties of periodic crystalline phases are predicted entirely 
from first principles. These approaches are especially at¬ 
tractive for predicting the properties of systems that are 
otherwise too expensive or difficult to study experimen¬ 
tally, such as alloys containing Rrt^, To^, and Pfi, among 
others^-—. 

While remarkably successful in predicting phase sta¬ 
bility in a wide variety of chemically disparate systems^, 
the occasional failures of DFTii^rii highlight the impor¬ 
tance of experiments to validate such predictions. How¬ 
ever, there exists a fundamental challenge in comparing 
DFT and experiments: electronic structure calculations 
predict zero kelvin properties, while experiments are per¬ 
formed at finite temperatures. The most reliable mea¬ 
surements of thermodynamic properties are performed at 
elevated temperatures, where equilibrium is more read¬ 
ily attained, but also where the entropic contributions to 
such properties are the greatest. Accuracy in compar¬ 
ing ab-initio and experimental results is vital, because 
any mismatch may indicate failure in the approxima¬ 
tions used in DFT to accurately reproduce the necessary 
physics. 

Of particular concern are disagreements between the 
set of observed phases and DFT-predicted zero kelvin 
ground states. These errors, while more sensitive to un¬ 
known kinetic barriers or unaccounted for entropic con¬ 
tributions, can also indicate fundamental flaws in the 
ab-initio method. Binaries that pair late 3d with late 
bd transition metals, such as Cu-A n^^i^^ , Co-Pt^iii, Ni- 
Fe-Ptiii^, and Fe-Pd^S, represent one class of 


materials where DFT predicts a rich variety of zero kelvin 
ground states for which no experimental evidence ex¬ 
ists. A wealth of long-period superstructures have been 
predicted to be stable in the intermediate continuum of 
compositions between xpt 0.5 and 0.75 in the Nii_a,Pta^, 
Fei_a;Pta;, Coi_a^Pta^ and Cui-a^Aua^ alloys, instead of the 
two phase mixtures of LIq and LI2 observed experimen¬ 
tally. In the case of Cu-Au, Co-Pt, Ni-Pt and Fe-Pd, the 
LI 2 AB 3 structure is altogether excluded from the set of 
ground states, and in the case of Co-Pt, the LIq forma¬ 
tion enthalpy is less stable than the (experimental) solid 
solution enthalpy21. In all cases, the formation enthalpies 
for the ordered phases have been predicted to be dozens 
of meV higher than experimental results^”—. These re¬ 
sults derive from zero kelvin predictions, but little or 
no thermodynamic analysis of the finite-temperature im¬ 
pacts has been performed, and bulk phase diagrams de¬ 
riving solely from electronic structure calculations have 
never previously been constructed for these materials. 
Co-Pt, Fe-Pt, and Fe-Pd are all candidates for use in 
ultrahigh density magnetic storage^”—; resolving un¬ 
certainties about low-temperature predictions of phase 
(in)stability is thus critical. 

Here, we explore the finite temperature implications 
of the zero kelvin ground states predicted by DFT for 
the Coi-x^^x alloy. We have developed a first-principles 
cluster-expansion Hamiltonian and used it in semi-grand 
canonical Monte Carlo simulations to construct a tem¬ 
perature versus composition phase diagram. The pre¬ 
dicted phase stability is in qualitative disagreement with 
available experimental observations. Our calculations 
also indicate that long-range ordered phases persist as 
ground states when including corrections for spin-orbit 
coupling, antiferromagnetic and mixed magnetic order¬ 
ing, and noncollinear magnetism. Using a hybrid func¬ 
tional, we have attempted to recover experimental re¬ 
sults; however, this approach is shown to introduce new 
errors in the enthalpy and magnetic moments of the 
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structures. These errors can be traced back to known 
failings of Hartree-Fock exchange when applied to tran¬ 
sition metals, which we illustrate by analyzing the density 
of states (DOS) for LIq CoPt. 

II. METHODS 

A. Electronic Structure Calculations 

All electronic structure calculations were performed 
using the Vienna Ab Initio Simulation Package 
(VASP)^“— with projector-augmented wave (PAW) 
potentials^i^ using the PBE functional^ for exchange 
and correlation (XC) energies^. In all cases, the en¬ 
ergy cutoff for the plane wave basis set was 460 eV and 
a F-centered Monkhorst-Pack^ A:-point mesh, converged 
to energy changes of less than 1 meV/atom, was used. 
All simulation cells were allowed to relax their volume, 
atomic positions, and collinear magnetic spins. Spin po¬ 
larization was employed and the correlation interpolation 
formula of Vosko et al^ was used to enhance magnetic 
energies. 

B. Cluster Expansion and Monte Carlo Simulations 

Configurational degrees of freedom in an alloy can be 
described with a cluster expansion Hamiltonia n^^d^ com¬ 
bined with Monte Carlo simulations^”—. The fully re¬ 
laxed energy [E(d)] of the arrangement [a] of two com¬ 
ponents on a parent crystal structure can be expressed 
via the following Hamiltonian)^ 

i;(a) = yo + V^« E w 

a Se^oc 

where a indexes a type of cluster (e.g., a nearest- 
neighbor pair, a next-nearest-neighbor pair, a nearest- 
neighbor tetrahedron, etc.), is the set of clusters 
symmetrically-equivalent to a, and 6 indexes different 
clusters in . (ps (d) is the cluster function for the cluster 
indexed by S, and the Va are effective cluster interactions 
(ECIs) obtained from projecting (e.g., via regression) ab 
initio energies and configurations onto the basis set of 
clusters (with Vb as the “empty cluster”, i.e., a constant 
term). We chose values for of -1 for Co and +1 for Pt, 
and 02 (d) = product of the of the 

sites in the cluster). 

Semi-grand canonical Monte Carlo simulations were 
performed in a 24 x 24 x 24 periodic supercell with 1000 
equilibration passes followed by 2000 thermodynamic¬ 
averaging passes (where one pass is NgUes attempted 
“spin flips”). Approximate phase boundaries were iden¬ 
tified by discontinuities in the composition-temperature 
lines at constant chemical potential (using temperature 
increments of AT = 2 K), or by plateaus in composition- 
chemical potential lines at constant temperature (using 
chemical potential increments of A/i = 0.01 eV). 


C. Phonons 

Eorce constants were calculated using the frozen 
phonon approach^”—, perturbing large super cells (108 
atoms) with small, isolated, symmetrically-distinct 
atomic displacements ( 0 . 01 , 0 . 02 , and 0.04 A). Eollowing 
electronic structure calculations, the resulting force con¬ 
stants were used to construct the dynamical matrix^. 
Vibrational free energies^ were then calculated using 
the quasi-harmonic approximation, repeating the previ¬ 
ous procedure at a variety of slightly smaller and larger 
(— 2 % to + 10 %) supercell volumes and using a second 
order polynomial to fit the dependence of the free energy 
on volume to determine the change in formation energy 
with temperature^. 


III. RESULTS AND DISCUSSION 

A. Zero kelvin and Finite Temperature Results for 

PBE 

To fully characterize the ab-initio properties of the 
Co-Pt binary, we calculated the DFT energies of 1469 
symmetrically-distinct orderings on the FCC lattice se¬ 
lected by an iterative approach. Each configuration 
was initialized ferromagnetically during the DET calcu¬ 
lations. We started with all known ECC based ground 
states as well as all symmetrically distinct orderings of 
Co and Pt over ECC within supercells containing up to 
6 atoms. A cluster expansion was fit to these energies 
and subsequently used to search for low energy config¬ 
urations in larger super cells. The energies of configura¬ 
tions in the larger supercells that were predicted to be 
below or close to the convex hull with the cluster expan¬ 
sion were then calculated with DET and included in a 
new fit. The cluster expansion was iteratively improved 
until no new ground states were predicted. The resulting 
set of structures included: (i) all unique supercells up 
to 6 atoms, (ii) all supercells up to 8 atoms with plat¬ 
inum composition between 25% and 75%, and (iii) long- 
period superstructure orderings involving ( 001 ) layers of 
Co and Pt stacked in various tilings in supercells with 
up to 15 atoms (1 x 1 x n primitive cells). Additional 
magnetic configurations were tested; anti-ferromagnetic 
(AEM) and mixed/ferrimagnetic orderings were explored 
(details of the orderings are given in the Supplemental 
Material, Eig. SI—), and non-collinear magnetism with 
spin-orbit coupling (SOC) corrections was separately in¬ 
cluded. The formation energies of all these configura¬ 
tions, calculated using HCP Co and ECC Pt as reference 
states, are shown in Eig. [H 

The set of ground states predicted by PBE (red 
circles/lines in Eig. m\ corresponding structures in 
Eig. |l(b)[ ) contain a variety of ordered phases at differ¬ 
ing compositions. Most prominent is the (32 CoPt 2 or¬ 
dering, characterized by alternating single (001) layers 
of Co and pairs of (001) layers of Pt. The LIq struc- 
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FIG. 1. (color online) (a) Formation energies in the Co-Pt bi¬ 
nary as calculated using PBE and HSE06. Eor ground state 
configurations, antiferromagnetic and mixed/ferrimagnetic 
orderings and SOC corrections were also considered. Exper¬ 
imental values at 914K are from Ref. [H- (b) Crystal struc¬ 
tures of the stable ground states predicted in Co-Pt using PBE 
(drawn with VESTA—). Compositions for ground states, in¬ 
cluding long-period superstructures, are indicated by dotted 
lines. 


ture, observed experimentally through a wide range of 
compositions^i^, is predicted to be stable only in a very 
narrow chemical potential range as a result of the small 
difference between the slopes of the tie lines connecting 
Llo CoPt to Co 2 Pt and Co 2 Pt 3 . The experimentally ob¬ 
served LI 2 CoPts ordering is entirely absent from the 
set of ground states, excluded by the depth of the com¬ 
mon tangent between the CoPt 2 and FCC Pt phases. A 
large number of long-period superstructures, character¬ 
ized by different arrangements of (001) layers of Co and 
Pt and with a preference for Pt-Pt layer pairs (similar to 
P 2 ) to account for deviations in stoichiometry (instead of 
anti-site defects on a LIq supercell), are present along or 
within 5 meV of the tangent between CoPt and CoPt 2 . 


These enthalpies and ground states match the results 
of Chepulskii et al^ and Hart et although an or¬ 
der of magnitude more configurations have been consid¬ 
ered here. While both prior works have reported a stable 
D 0 i 9 (HCP) phase at Co-rich compositionsiiii, we re¬ 
stricted our focus to FCC-like superstructures because 
we are only interested in equiatomic and majority-Pt al¬ 
loys (i.e., 0.25 < xpt < 0.75), which all adopt FCC-based 
orderings experimental!}^. 

The calculated formation energies of LIq CoPt and LI 2 
CoPts are approximately 50 meV above (less stable) the 
measured formation enthalpies at 914 K— , i.e., an er¬ 
ror of nearly 40%. The lattice parameters, however, are 
within <1% of experimental measurements^. Alterna¬ 
tive magnetic configurations and spin-orbit contributions 
did not result in any further lowering of the formation en¬ 
ergies of the ground states and increased the difference 
between calculated and measured formation enthalpies 
[blu e tria ngles, green squares, and purple pentagons in 
Fig.[U a)| respectively]. 

The Co-Pt system forms a FCC-based solid solution at 
high temperatures and a variety of FCC-based ordered 
structures upon cooling. The contributions of configu¬ 
rational entropy play an important role in determining 
phase stability when increasing temperature. We fit^ 
the coefficients of a cluster expansion to the 1469 for¬ 
mation energies using a genetic algorithm^ to determine 
the optimal basis set. The resulting fit has 89 effective 
cluster interaction (ECIs) coefficients corre spond ing to 
pair, triplet, and quadruplet clusters [Fig. |2(a)| . The 
root-mean-squared error of the fit was 3.3 meV, and the 
cross validation score (using leave-one-out cross valida¬ 
tion) was 3.6 meV. 

The PBE-based phase diagram [Fig. |2(b)] , determined 
using semi-grand canonical Monte Carlo, shows a very 
wide stability region for (32 CoPt 2 . Although this order¬ 
ing also has the highest order-disorder temperature, the 
transformation is predicted to occur hundreds of degrees 
lower than the experimental transition temperatures of 
Llo CoPt and LI 2 CoPts— The LI 2 CoPts struc¬ 
ture does not appear at all in the calculated phase dia¬ 
gram, and the region surrounding xpt = 0.5 consists of 
a continuum of defected incommensurate long-period su¬ 
perstructures up to the peritectoid temperature, decom¬ 
posing into a mixture of solid solution and P 2 CoPt 2 . 
We found no evidence for the stabilization of defected 
Llo CoPt (i.e., single planes of Co and Pt with a random 
distribution of anti-site defects); even at compositions 
close to xpt = 0.5, the structure instead resembled /32-like 
CosPts with anti-site defects concentrated in the paired 
Pt-Pt (001) layers. The phase diagram was not explored 
below 400 K, nor were compositions below xpt = 0.2 or 
above xpt = 0.8 explored because these regions were not 
emphasized when fitting the cluster expansion. 

Overall, the resulting phase diagram is inconsistent 
with the high temperature observations of LIq CoPt and 
LI 2 CoPts. These phases have been well characterized in 
the literature and have important differences in diffrac- 
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FIG. 2. (color online) (a) Effective cluster interactions of 
the PBE cluster-expanded Hamiltonian, grouped by cluster 
type. The empty and point clusters have ECIs of 20.5 meV 
and 473 meV, respectively, (b) The temperature-composition 
phase diagrams for Co-Pt from experiments (dashed red, blue 
lines^i^i^), our results with PBE (green), and CALPHAD 
(solid black lines^). Green shading indicates two-phase re¬ 
gions. 


tion patterns from the long-period superstructures pre¬ 
dicted by PBE. Experimental characterization of the LIq 
phase has historically relied upon the [001] superstruc¬ 
ture peab^i^i^, a peak that is absent in the predicted 
long-period superstructure orderings. LI 2 CoPts, having 
cubic symmetry, also has very different diffraction pat¬ 
terns from tetragonal (32 CoPt 2 or any of its derivatives. 

While the calculated phase diagram only accounts for 
configurational excitations at finite temperature, vibra¬ 
tional and magnetic excitations may also play a role 
in determining the relative stability of different order¬ 
ings. To explore vibrational free energies, we performed 
phonon calculations for FCC Co, FCC Pt, LIq CoPt, 
(32 CoPt 2 and LI 2 CoPts within the quasi-harmonic ap¬ 
proximation. However, even at 914 K, LI 2 CoPts never 
emerges as a stable phase relative to /32 CoPt 2 and FCC 
Pt. Furthermore, the stability of LIq CoPt relative to (32 
CoPt 2 and LI 2 CosPt does not increase markedly with 
increasing temperature. This suggests that rigorous in¬ 
clusion of vibrational degrees of freedom, together with 


configurational degrees of freedom using coarse-graining 
schemes^i^, are unlikely to qualitatively alter the cal¬ 
culated phase diagram of Fig. |2(b)[ 

Thermal excitations of magnetic moments have previ¬ 
ously been shown to be important in affecting phase sta¬ 
bility in Co-Pt alloys^^A^. In pure Co, entropic contribu¬ 
tions arising from spin-spiral excitations have a strong in¬ 
fluence on the HCP/FCC transformation temperature^. 
At Pt-rich compositions, the ferromagnetic to para¬ 
magnetic Curie temperatures occur below the order- 
disorder transition temperatures of both LIq CoPt and 
LI 2 CoPtgrliSi^. Hence, magnetic disorder, neglected 
in the calculation of our phase diagram, will likely play 
a role in determining the precise order-disorder transi¬ 
tion temperatures. An in-depth study of the effects of 
magnetic entropy on the phase diagram are beyond the 
scope of this study. However, below the Curie temper¬ 
atures contributions to the free energy from magnetic 
entropy are expected to be small. While we do not know 
the Curie temperature of the (32 phase, those of LIq and 
LI 2 are above all of our calculated solid solution transfor¬ 
mation temperatures. In the case of the (32 phase, there 
are two possibilities: (i) CoPt 2 remains ferromagnetically 
ordered up to the solid solution transformation tempera¬ 
ture, or (ii) CoPt 2 becomes paramagnetic below the solid 
solution transformation temperature. In scenario (i) our 
phase diagram should be negligibly impacted by the in¬ 
clusion of magnetic entropy, while in scenario (ii), mag¬ 
netic entropy will only further stabilize the ^2 phase with 
respect to the LIq and LI 2 phases. We therefore expect 
that inclusion of magnetic excitations will not rectify the 
disagreement between finite temperature predictions and 
experimental observations. 

The over-stabilization of (32 and other long-period su¬ 
perstructures with paired (001) Pt-Pt planes can be at¬ 
tributed to over-delocalized electron charge densities in 
PBE. Generalized gradient approximations (GGAs) can 
perform poorly in transition metals^^— where the sig¬ 
nificance of the localized d-orbitals comes into conflict 
with the orbital-less approach of GGA. Any functional 
treatment of the electron density leads to an electron 
interacting with its own potential, the self-interaction 
error, which causes excessive delocalization of the total 
charge. Additionally, the GGA functional can not ener¬ 
getically differentiate between occupied and unoccupied 
bands, leading to incorrect predictions for orbital/band 
occupation and splittin g^^i^^ . Because both the energy 
levels and occupations of the d orbitals are incorrect, the 
Co 3d and Pt 5d bands cannot hybridize, losing signifi¬ 
cant enhancement of the magnetic moment^-—. Stabi¬ 
lization of a magnetic ground state, however, is a driving 
force in choosing the thermodynamic ground states in 
Co-Pt^, and ferromagnetic effects drive the asymmetry 
in the phase diagram with respect to the Co-rich and 
Pt-rich LI 2 phases^i^. 
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B. Comparison with Hybrid Functional HSE06 

Using a screened form of Hartree-Fock exchange, 
Zhang et alr^ were able to recover the experimental 
ground-states of Au-Cu. The parametrization of the hy¬ 
brid functional developed by Heyd et al. (HSE06)^i^ 
reduces the self-interaction error akin to Hartree-Fock 
theory^, while avoiding the associated singularity in oc¬ 
cupation at the Fermi level. Since HSE06 explicitly in¬ 
cludes orbitals, d-orbital hybridization can be recovered, 
e.g., as in Au-Cu. The functional is, however, limited by 
its computational expense^ and accurac y^^i^^ . 

To examine the performance of the HSE06 functional, 
all of the PBE ground states, as well as LI 2 CoPts, 
FCC Co and Pt, and HCP Co were recalculated [pink 
hexagons/lines in Fig. |l(a)| . The results show qualita¬ 
tive improvement: LI 2 CoPts is predicted as a ground 
state, and LIq CoPt is substantially more stable relative 
to LI 2 CosPt and ^2 CoPt 2 . However, the enthalpy and 
lattice parameter errors increase by an order of magni¬ 
tude (200-300 meV and 2.5-4.5%, respectively) and ^2 
remains as a ground state. These results are not surpris¬ 
ing: HSE06 is known to severely overestimate exchange 
splitting in itinerant magnetic systems2^“— . The default 
screening parameter in the HSE06 formalism {uj = 0.2 
screen = 10 A) results iu an effective screening 
length an order of magnitude greater than the length 
scales {rscreen = 0.24 — 0.26 A) for screening in bulk 
Pt or Co metals^^, introducing spurious interactions be¬ 
tween orbitals at different sites. This is further verified 
by stabilization of a ferromagnetic ground state for Pt in 
both this work and Ref. over the experimentally ob¬ 
served nonmagnetic state, and an overestimation of mag¬ 
netic moments with HSE06 for all experimentally studied 
structures. 

The shortcomings of both the PBE and HSE06 func¬ 
tionals are more easily visualized by comparing the cal¬ 
culated density of states (DOS) for Co-Pt (Fig. [3j) with 
data from experiments. The width of d-band states, 
as measured by photoemission^UiSi, are shown in grey 
in Fig. O while the electron occupation at the Fermi 
level can be inferred from low temperature heat capac¬ 
ity experiments^-—, shown as the height between the 
red dashed lines. In the case of HCP cobalt, PBE pre¬ 
dicts a d-band width in good agreement with experiment, 
but the occupation at the Fermi level is underestimated. 
The excessive delocalization provided by self-interaction 
errors in PBE smears the electron density towards a more 
even distribution between orbitals, and this effect is mag¬ 
nified by the inability of PBE to (energetically) distin¬ 
guish between occupied and empty orbitals. This latter 
point is exemplified by a difference between the calcu¬ 
lated exchange splitting (energy difference between max¬ 
ima in spin-up and spin-down DOS) and the experimen¬ 
tally determined width, indicated in purple. Unsurpris¬ 
ingly, PBE performs well in FCC platinum, where the 
d-orbitals are fully occupied and differentiation between 
sub-shells does not matter. Though we were not able 


to find suitable experimental analysis of the LIq CoPt 
d-bands. Fig. [3] shows an unexpectedly wide d-band and 
low electronic occupation at the Fermi level in CoPt, sim¬ 
ilar to Co. 


Exp. Exch. 



FIG. 3. (color online) d-orbital density of states for HCP Co 
(top), FCC Pt (middle) and Llo CoPt (bottom) calculated 
from DFT, using both the PBE (blue line) and HSE06 (green 
line) functionals. Gray regions show occupied d-band width 
from photoemission results^i^, red dotted lines show occupa¬ 
tion at the Fermi level predicted from the heat capacit}^"— , 
and purple dotted lines show exchange splitting predicted 
from experiment [AE between for majority and mi¬ 

nority spins] for Co^. 

Contrary to expectation, HSE06 worsens the degree of 
electronic delocalization in Co even as the self-interaction 
error should be reduced; this can be attributed to the 
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unrealistic (for Co) screening length standardized in the 
HSE06 functional. In Pt, a very large fictitious exchange 
splitting is introduced, explaining the large magnetic mo¬ 
ment commented upon previously. Here, an electron has 
been promoted from the bd orbitals in PBE to the 6p or¬ 
bitals in HSE06 to increase the net number of unpaired 
spins. The over-stabilization of a magnetic ground state 
also reduces the energy of the d-orbitals in all materials; 
we can see this in Eig. [3] where the occupation at the 
Eermi level drops to between one half and one quarter of 
the PBE values. These results both highlight the failures 
of the PBE and HSE06 functionals and offer insight as to 
why these methods describe the Co-Pt binary so poorly. 


C. Additional Functionals 

Additional functionals were explored, though not in- 
depth, once initial results indicated a qualitatively simi¬ 
lar set and arrangement of ground states to the results us¬ 
ing PBE. The Tao-Perdew-Staroverov-Scuseriai^ meta- 
GGA functional and its revised variant^, as well as the 
Minnesota meta-GGA functional MOG-L^, were tested, 
but retained / 32 -like structures and excluded LI 2 CoPts 
from the set of ground states. PBESol^ was investigated 
but also produced similar results to PBE. HSESol2£, and 
various parametrizations of the HSE functional (in a and 
uj) were tested only on LIq GoPt, LI 2 G 03 Pt and GoPts, 
and (32 GoPt 2 , but yielded decreasingly tiny occupations 
at the Eermi level as well as worsening energetics and a 
deepening of the (32 enthalpy with respect to other ground 
states. 

The importance of localized d electrons in Go-Pt, 
which has been incorrectly treated in both PBE and 
HSE06, motivated the examination of the LSDA+U 
methodSl. We explored a two-dimensional grid of I/’s 
of 1.0-4.0 eV in Go, and 0.0-4.0 eV in Pt, in 0.1 eV 
steps using the rotationally-invariant method of Dudarev 
et al^. Unfortunately (though perhaps not unsurpris¬ 
ingly), the results showed no set of I/’s that simultane¬ 
ously produced the correct set of ground states, matched 
the experimental enthalpies (within a range of ±50%), 
and matched the experimental c/a ratio (within a range 
of ±50%). The linear-response approach to determine 
coupling parameters of Cococcioni and Gironcoli^^ was 
also used, but the resulting I/’s produced similar unsat¬ 


isfactory results. 

IV. CONCLUSION 

In summary, we have shown that the anomalous 
ground states predicted in Co-Pt using the PBE func¬ 
tional result in equally anomalous phase behavior at ele¬ 
vated temperature, wholly inconsistent with the body of 
experimental literature. By examining our results in the 
greater context of the known shortcomings of PBE, we 
can characterize the modes of failure, attributing the sta¬ 
bilization of the (32 CoPt 2 ground state to self-interaction 
and occupation errors inherent in the functional. Al¬ 
though DET performs exceptionally well in a wide vari¬ 
ety of inorganic systems, caution must be used when pre¬ 
dictions appear inconsistent with experiment. By using 
rigorous statistical mechanical approaches, experimental 
results can be meaningfully compared with zero kelvin 
predictions in both qualitative and quantitative fashions. 
Unfortunately, when one method falls short, it is not al¬ 
ways sufficient to move up the “Jacob’s ladder” of XC 
functionals^^: the HSE06 functional merely trades one 
set of inaccuracies for another. Based on both our analy¬ 
sis and the existing literature, we believe that sibling sys¬ 
tems (e.g., Ee-Pt, Ni-Pt, Ee-Ni) will yield similar results. 
Though the prospect of a one-size-fits-all DET-based ap¬ 
proach to predicting phase diagrams is appealing, our 
analysis of the Go-Pt system highlights the need to re¬ 
view zero kelvin electronic results in a finite-temperature 
and thermodynamically meaningful fashion. 
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